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Abstract 

We generalize Hamilton's principle with fractional derivatives in La- 
grangian L{t,y{t),oD"y(t),a) so that the function y and the order of 
fractional derivative a are varied in the minimization procedure. We de- 
rive stationarity conditions and discuss them through several examples. 
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1 Introduction 

Hamilton's principle is one of the basic principles of Physics. Anthony [4] states: 
"In theoretical physics, a theory is often considered to be complete if its varia- 
tional principle in the sense of Hamilton is known" . When a Hamilton's principle 
is known, the whole information concerning the processes of a particular system 
is included into its Lagrangian. When known, the Hamilton's principe could be 
used in many different ways. For example, it could serve as the basis for ob- 
taining first integrals via Nother's theory, or to generate approximate solutions 
to the relevant system of equations by the use of Ritz procedure (cf. [31]). 

In this paper we investigate necessary conditions for solutions of fractional 
variational problems (Euler-Lagrange equations). Such investigations has been 
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initiated in [351 121]) and continued in [T] - [S] and [5]; see also [21] for the 
importance of introducing fractional derivatives into the Lagrangian density of 
Hamilton's principle. In general, we refer to [5] - [IB ] . [TOll^l^ . [^ - [5^. for 
different aspects of the calculus of variations and fractional calculus, motivations 
and applications. 

When a fractional variational problem is studied, a natural question arises 
how one can choose a, the order of fractional derivative, in order to achieve the 
minimal value of a functional under consideration. More generally, one can ad- 
dress the same question for any problem involving fractional operators. Usually, 
in application, a good choice of a is obtained by experiments, numerical meth- 
ods or computational simulations. However, experimental results give different 
values for a within certain interval. In this paper we propose a method which 
a priori gives values for a which optimize the considered variational problem 
following the fundamental minimization principle of Hamilton's action. In fact, 
we address the question of finding stationary points for Hamilton's action in- 
tegral with fractional Lagrangian in a more general setting. Namely, we allow 
the stationarity of the action integral with respect to a set of admissible func- 
tions and with respect to the order of fractional derivatives, appearing in the 
Lagrangian. Up to our knowledge the problem when both y and a are varied 
has not been analyzed so far. It leads to stationarity conditions as a basis for 
generalized Hamilton's principle for the action integral 

I[y,a]~[ L{t,y{t),oD^y,a)dt, y Gl(,a G A := [O,ao],ao < 1, (1) 
Jo 

where W is a set of admissible functions: Find 

min ^[y, a] or (2) 

{y,a)eUxA 

min(min/[y, q;]) or (3) 

min(min/[?/, a]). (4) 

In this paper we analyze stationarity conditions for ([2]) and ([3]). Stationarity 
conditions with respect to a in ([4]) are more difficult, and, contrary to (|3|), this 
case is less natural in applications. Note that in ([5]), (O and (|3]) one can look 
for maximums instead of minimums. So, our general problem is determination 
of stationary points. 

So far, parameter a, the order of fractional derivative, has been determined 
experimentally (cf. e.g. [37]). This approach offers a rational way for choosing 
the precise a. 

The paper is organized as follows. To the end of Introduction, we recall 
definitions and properties of fractional derivatives. In Section [2] we present a 
framework in which we shall study variational problems ([2]) and (|3]). Then in 
Section [3| we derive stationarity conditions for ([2]). Section |4| is devoted to addi- 
tional assumptions which provide equivalence of problems ([2]) and (|3]). Results 
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which are obtained in previous sections are illustrated by several examples in 
Section [5] Moreover, examples of this section give further motivation for our 
investigation. In the last remark of Section [S] we propose a new formulation of 
a fractional variational problem. 

Throughout this paper we shall use the following notation. The mapping 
{t,a) i-> o-Df (?/), which defines the left Riemann-Liouville fractional derivative 
of order a, will be denoted by oD^y, or by aDfy{t). Recall, 

where V is the Euler gamma function, and its existence is provided whenever 

m^t^jl^^dr (5) 

is an absolutely continuous function. Recall, the space of absolutely continuous 
functions is denoted by ylC([0, 6]) and it is supplied with the norm ||/|| = 
sup3.g[o,6] 1/(2^)1 (clearly, it is not a Banach space). For example, ([5]) is absolutely 
continuous if y G AC([0, 6]). However, there are some cases when with less 
regularity in y we still have a well-defined operator of fractional differentiation 
(cf. [2l])- For instance, oDfy exists for functions with integrable singularities (a 
continuous and locally integrable function / in (0, 6] has an integrable singularity 
at T = of order r < 1 if lim^ _>o T''/(r) ^ 0). In particular we can take 
y{t) = i"'', t G (0, h] (for any 6>0),0<^<1. Then we obtain the so-called 
Euler formula (cf. ^ (2.26)]) 

oA"^'" = J^}^ \ , te(0,6]. 

The right Riemann-Liouville fractional derivative of order a is defined as 

The conditions for its existence are similar as in the case of the left fractional 
derivative. 

In the sequel we shall consider cases involving both fractional derivatives 
and work with integrable functions for which these derivatives (or one of them) 
exist. In such cases notation oD^y, resp. tD^y, t G [0,6], means that y and 
QD"y, resp. tD^y, are considered as integrable functions which can take values 
-foo or — oo at some points. 

Recall (gl]), that oDfy y' and jD^y ~y' in C([0,6]), as a ^ 1", 
whenever y G C^([0, h]). 

Also, we shall make use of Caputo fractional derivatives. The left, resp. 
right, Caputo fractional derivative of order a G [0, 1) is defined as 

CDty:=—^,f-^^^dT resp. f := _1_ /' dr. 

° r(l-a)io (i-^)" t by T{l-a)Jt (r - 1)" 
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One can show that for y E AC{[0, b]) and t G [0, b], 



oD'^y = ^D^y + — TT-^^, tD^y^'i'D^y 



r(l-a) (t-a)"' "'^ * r(l-a)(6-i)"' 

(cf. e.g. [12). Therefore, o-D"y — ^Dfy, resp. t-D^y = ^D^y, whenever 
2/(0) = 0, resp. y(&) = 0. 

2 Formulation of the problem 

We investigate stationary points of ([T]) for a G [0, ao] and all admissible functions 
y, whose properties will be specified in the sequel. We shall distinguish two cases: 
ao strictly less than 1 and ao = 1. In the case ao < 1, set 

Ui := {y e L\[0,b])\oD^y e L\[0,b])}. 

Obviously, AC{[0, b]) is a subset oiUi. In the case ao = 1 we assume that y ElAi 
and that, in addition, oD^y exists, and oD^y = y' is an integrable function. Let 
us note that one can consider Ui defined with LP{[0,b]) (or their subspaces) 
instead of i^([0, b]) (see Remark [23]). 
In general, we shall use notation 

U :~ {y E Ui \ y satisfies specified boundary conditions}. (6) 

We shall sometimes write U also for Ui (then the set of specified boundary 
conditions is empty). 

In the sequel, Lagrangian L{t, y(t), ^Dfy^ a) (Lagrangian density, in Physics) 

ieCi([0,6] xRxRx [0, 1]) \ 

and \ (7) 

t K> d:iL{t, y{t), oD^y, a) eUr, for every y eUi J 

where Ur {y e iH[0,&]) I tD^y £ L\[0,b])}. 

Recall, our generalization of Hamilton's principle is realized through the 
determination of (t/*,a*) €U x A such that 

min / L{t,y{t),oD^y,a)dt= / L{t,y*{t), oOf y* , a*) dt. (8) 

There are two special cases of ©. The first one is obtained when A = {1}. 
Then, since o-Dfylo^i = y'{t) for y £ C^{[0,b]), the solution y* of dH) satisfies 
the classical Euler-Lagrange equation 

d dL 9L 
dt dy' dy 

If A has a single element A = {a},0<a<l, then min I[y,ct] leads 

(y,a)eUx{a} 

to the fractional Euler-Lagrange equation (cf. ([11 IS])) 
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We proceed with finding stationary points related to ([T]). 



3 Optimality conditions 

A necessary condition for the existence of solutions to variational problem ([8]) 
is given in the following proposition. 

Proposition 3.1 Let L satisfy Then a necessary condition that functional 
(OP has an extremal point at {y*,a*) ^ U x A is that (y*,a*) is a solution of 
the system of equations 

f+^D^^ = 0, (9) 
— Giy,a) + —]dt = 0, (10) 



\doDfy ' da 
where 

G{y,a) = = ±(f,^^y)(t,a), f,{t,a) = -^J—--[^{l-a)-lnt], t > 0, 

with the Euler function ip{z) = Inr(z), and (/i *t y)it, a) — J* /i(t, a)y{t — 
t) dr. 

Proof. Let {y* , a*) be an element oi U x A for which I[y, a] has an extremal 
value. Let y{t) = y*{t) + eif{t), a = a* + £2, £1,62 G M, with / G Ui, and 
the boundary conditions on / are specified so that the varied path y* + Sif is 
an element of ^. Then I[y,a] — I[y* +£1/, a* +£2] =: /(£i,£2)- A necessary 
condition for an extremal value of I[y, a] is 



a/(£i,£2) 



dei 

Therefore we obtain 



Q 9/(£l,£2) 



= 0. 

£1=0, £2=0 



Applying the fractional integration by parts formula (cf. [17]): 

g{t) ■ oA"/(i) dt = / fit) ■ tD^g{t) dt, (13) 



/o ^0 
to ((TT|) . (fTTj) is transformed to 
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From this equation, using the fundamental lemma of the calculus of variations 
(see [9l p. 115]), we conclude that condition ^ holds for the optimal values y* 
and a* . In (jlOp the term ^ is transformed by the use of expression 



da ' ' T{l-a)dtJo {t ~ t) 

= j^{h*ty){t,a) 

= Giy.a), {y,a)eUxA (14) 
(cf. O p. 592]). We obtain ^ by substituting 1^ into 1^^. □ 

Remark 3.2 In general, in solving equations ^ and (fTO|). the most delicate 
task is the calculation of expression ^ . Although its general form (|14l) has 
been derived in [SI p. 592], various difficulties can still appear. We illustrate 
this by examples in Section [SJ However, the simplified form of ^ , in some 
special cases, is important. 

Already in 6 it has been shown that for y e AC([0, b]) 

= -{-f + \nt)y{0)~f\-f + lnT)y{t-T)dT 



da 



Q=0 + 



i, + lnt)yit)+ f y^'^- y^'-^^ dr, (15) 



IS 

Q=0 + 



' 

where 7 = 0.5772156... is the Euler constant. (Another form of ^"q^ ^ 
also given in [32j p. 111].) 

Let us obtain a simplified form of ^°g^ ^ at a = 1^. In order to do that we 
use the method proposed in [53]. We recall the expansion of {t — tY /r{l + e) 
with respect to e, at e = 0, with t < t (cf. [IHl p. 401]), which will be used in 
the sequel: 



= l + e(7 + In(i-r))+o(e). (16) 



r(l+e) r(l + £) 
Assume now that y G C^([0, 6]). Then, as in [23j, for t e [0,6], 



r(l - a)i" r(l - a) Jo {t - t) 

y(0) , y^'Ho) , 1 /•* 



r(l - a)i" r(2 - a)t«-i r(2 - a) Jo 
With a = 1 — e we obtain 
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From ((ni) and (HU) it follows that 



da 



doD}- 



t 



=0+ 



y(i)(0)lnt-72;«(<)- 



2/(2) (^) ln(t - t) dr. (18) 



Assuming 2/(0) = in (fT8|) . we recover the resuhs presented in 29 and [30 for 
the Caputo fractional derivative. Since oD"y = oD^y when y{0) = 0, it follows 
that with 2/(0) = ([18]) becomes 



doDfy 



da 



d'^Dfy 



da 



-2/(^HO)lni-72/(^'(i)- / 2/^^' M In(t-r) dr. 



Remark 3.3 Functional ([T|) is a special case of the class of functionals with La- 
grangians depending on linear operators, see |12[ p. 51]. Indeed, suppose that 
Lagrangian L depends on t, y and Cy, where C : M — i> i''([0, 6]), p G [1, +oo), is 
a linear operator acting on a set of admissible functions M, which is linear, open 
and dense in i^([0, h]) (i.e. C belongs to Lin{M, i^([0, 6])), the space of contin- 
uous, linear functions with the uniform norm). Suppose that L is continuously 
differentiable with respect to t and y and twice continuously differentiable with 
respect to Cy. Moreover, assume that function t i— L{t,y{t), Cy{t)), t G [0,b], 
is continuous, for all y G M. Then the Euler-Lagrange equation reads 



dL 
dy 



dL 



= 0, 



where C* denotes the adjoint operator of C In case when C is the left Riemann- 
Liouville operator oDf, with the adjoint fD^^ the latter equation coincides with 
®. 

If instead of C one considers a family {Ca,a G A}, where A = [0,1] (or 
some other interval), and the mapping A Lin{M, LP{[0,b])), a t—^ £„, is 
differentiable, then a more general problem of finding stationary points with 
respect to y and a can be formulated. In that case, one can derive the second 
stationarity condition similar to ()10|) : 



dLdC 
dC da 



dL 
da 



dt = 0. 



4 Equivalent problems 

In this section we shall give conditions which provide that problems ([2]) and p|) 
coincide. 

Proposition 4.1 Let the Lagrangian L satisfy Assume that for every a G 
[0,1] there is a unique y*{t,a) € U, a solution to and that the mapping 
a y*{t,a) is differentiable as a mapping from [0, 1] to U. Then the problem 
mi-a(^y a)£UxAl[y,o:] is equivalent to the problem minQGA(minye;^ /[?/, a]). 



7 



Proof. As we have shown in Proposition 13. 1[ any solution to the problem 
mini^y^a)eUxAl[y,ce] satisfies system ([9])- (fT0|) . It can be solved as follows. We 
first solve ^ and the corresponding boundary conditions to obtain y* = y* (t, a). 
According to the assumption, the solution y* is unique. Then we insert y* in 
PH]) to obtain a*. In this case, functional /[2/,Q;] becomes a functional depend- 
ing only on a, a 1-^ I[y*{t,a),a] = I[a], and therefore dTU]) transforms to the 
total derivative of I[a] since 



dl[a] 
da 



dl[a + e] 



de 



e=0 



dL 



dL dy 
dy da d^Dfy 



d_ 

da' 



dy 

da 



dL 
dy 



dL 



dL 



doDfy \da 



—oL>t 



doDfy 

y + 



dL 



d 



d^Dfy \da 



d£ 

da 
dL 
da 



dL\ , 
1 dt, 

da 



dt 



dt 



(19) 



where we used fractional integration by parts formula (|13p in the third, and 
equation © in the last equality. This proves the claim. □ 

The following simple assertion is of particular interest: 

Proposition 4.2 Let L satisfy Assume that for every a G [0, 1] there 

exists a unique G U, a solution to the fractional variational problem and 
that /[?/c(,a] is the corresponding minimal value of the functional /. Assume 
additionally that 

dl 



da 



Then the minimal, resp. maximal value of the functional I[y, a] is attained at 
a = 0, resp. at a = 1. 

Proof. Under the above assumptions we have that 

I[yo,0] < /[2/a,0] < I[ya,a] < I[yi,a] < Va e [0,1], 

which proves the claim. □ 



Remark 4.3 The same argument can be applied to the case when dl /da < 0, 
for any fixed ya G i.e. when / is an decreasing function of a, for any fixed 
ya gU. In that case the minimal, resp. maximal value of / is at a = 1, resp. 
a = 0. 
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5 Examples 

5.1 Examples with Lagrangians linear in y 

Example 5.1 Consider the action integral for the inertial motion (no force 
acting) of a material point of the form 



J[y,a]:- / {oDfyfdt, {y,a) x A, 
Jo 



(20) 



where U :={y GUi\ y{0) = 0, 2/(1) = 1} and A = [0, 1]. 

Obviously, the minimal value of I[y, a] is zero, and it is attained whenever 
o-Df 2/ ~ 0. Solutions to equation oD^y = arc of the form y{t) — C ■ t^^", 
t G [0, 1], C G R (cf. [H])- All solutions satisfy the Euler-Lagrange equation 

tD1{„D^y)^0. 

Stationarity condition PU)) reads: 

and is automatically satisfied. 

Note that C ■ t^"" G Ui, for aU C G M, but only t^-"' G U. Hence, we 
conclude that {y*,a*) — {t^~",a), a G [0,1], are solutions to the variational 
problem I[y,a\ — > min, for / defined by (I20|) . 

Remark 5.2 If L{t,y{t), oD^y, a) — {oDfyY + (a — ao)^, for a fixed ao G 
(0, 1), then the problem Lit, 2/(0: oDfy, a) dt — s> min has a unique minimizer 
(j/*,a*) = (ti-"«,ao). 

Example 5.3 Let the Lagrangian L be of the form 

L{t, y{t), oD^y, a) := {oDfyf ~c-y, c G E, 
and let 14 = {y eUi \ y{0) = 0}, A = [0, 1], for the variational problem: 

, ™in , liV: "] ^ , , / iioDfyf - c ■ y{t)) dt. 

Equations (|9]) and become: 



tD^D^y = c, (21) 

Equation (|211) could be solved as follows. First, one introduces a substitution 
2:(t) = oD'^y, t G [0, 1], and solve tDfz = c: 

z{t) = c- iG[0,l],aGA 

r(i - a) 
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Therefore, 



,D^y{t) = c-^^-^ te[0,l],a&A. (22) 

r(i - a) 

Recall, 

oI^y:=-^f {t~Try{T)dT, te[0,l],aeA 

r(a) Jo 

and apply it on the both sides of ([22]). Using Th. 2.4], i.e. oI^{oDfy) = y{t), 
one obtains 

t 



r(a)r(i + a) Jq 
- r(i + a)^4^^r(-a)p!r(i + p + a)^ ' ^^O'l]'"^^- 

This solution is unique and belongs to U. Since a H> a) is differentiable. 
Proposition 14.11 holds. 

We substitute obtained y{t, a) into I[y, a] which yields 



I\a] = 



Jo \\ r(l-a); ria)T{l + a)Jo ^ > ^ > j 
l-r^^2 c2 ^ r(p-a)r(l+p) \ 



r(l-a);' r(l + a)^^r(-a)p!r(l+p + a) 

Simple numerical calculations show that I[a\ is an increasing function and at- 
tains extremal values at the boundaries. 

Remark 5.4 Equation (PT|) represents a fractional generalization of the equa- 
tion of motion for a material point (with unit mass) under the action of constant 
force equal to c. Our result shows that an optimal value of Hamilton's action is 
attained for a = 1, that is for integer order dynamics. We note that different 
generalizations of classical equation of motion can be found in [TS] , where the 
problem oDfy = c, 1 < a < 2, was analyzed. 



5.2 Examples with Lagrangians linear in oDfy 
Example 5.5 Let 

L{t,y{t),aD'^y,a):^T{l-a)oD'^y-]^cy\ oQ.c^l, (23) 

and consider the problem of finding stationary points for the functional ([1]), 
where U := {y € Ui | y(0) = i} and ao < \- Note that L satisfies the so 
called primary constraint in Dirac's classification of systems with constraints 
(cf. [H]). In the setting of fractional derivatives such Lagrangians have been 
recently treated in 7 and [20j . 



10 



Equations ^ and ([TU]) become 

r{l - a)tD'^l - cy 



and 



' ,doD?y , ar(l-a) 

r(i - a,' 



dt = 0. 



(24) 



(25) 



da da 

Equation has a unique solution = di-t)" ' ^ ^ o; G [0, ao]- This 



implies 



1 



dt c(1 - T)°(t - r)" 2c(l-i) 
1 -1 



2a 



c(l-r)"(t-T)^ 



dr 







1 



1 



(l-r)2" 7o 2c(l-t)2" 



2c(l-t)2" 
dt 



dt 



1 



1 



2c Jo [i-ty^ 



dt. 



Since < we have that /[?/*, a] exists and is an increasing lunction with 
respect to a. Hence, I[y*,a] attains its minimal value at a = 0, and it equals 
l/(2c). We also have that the maximal value of I[y*,a] is attained at ao- 

Example 5.6 Let U := {y eUi\ y{0) =0}, cy^O and let L be of the form 

Lit,yit),oD?y,a):^c-oDfy + fiyit)), tG[0,l], (26) 

where the properties of / are going to be specified. 

In this example we are dealing with integrable functions which can take 
values +00 or — oo at some points. We are going to analyze stationary points of 



{c-oDfy{t) + f{y{t)))dt. 



Equations ([9]) and ([T0| become 



dl 

dy 



1 doDfy 



da 



dt = 0. 



(27) 
(28) 



Since tD^c = Y-{i-a){i-t)'' ' ^ ^ [OjI]' see that in order to solve ((27| - 
(PS)) we have to assume that / e C^(]R), and that /' is invertible so that 
i ^ (/0~^( r(i-aKi-t)° ^ ^ '^^^^ equation (P7|) is solvable with respect 
to y: 



yc{t,a) 



dl 
dy 



r(i -a)(i -ty 



t G [0,1]. 



(29) 
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Since c 7^ 0, dM]) implies dt = 0. Thus 

= f ^^^dt= f G{y,a){t)dt= C Ufi*tv){t,a)dt 
Jo oa Jo Jo dt 

= h*ty{t,a)\t=i- fi*ty{t,a)\t=o^ h*ty{t,a)\t=i, (30) 

where we have used that /i e i^([0, 1]) and that y £ U. Substitution of (^91) 
into (1201) gives (/i * yc{t,a)){t)\t=i = or 

/■^ V.(l-a)-ln(l-r) c \ _ n .o^^ 

io r(l-a)(l-r)" [dyj [t{1 - a){r - l)^ ) ^''^ 

Solving this equation is obviously difficult. Hence, we consider some special 
cases. 

a) f{y{t)) := d ■ t e [0, 1], deR. Then the Lagrangian is 

L{t, y{t),oDty, a)=c- oD^y{t) + d-^, 

and 

1 c 

yJt,a) = — -— — . ie[0, 11. 

' dr(l-a)(l-i)"' ^ ' ^ 

Also, ([3T|) becomes 

^i^(l_„)-ln(l-T) ^^^^ 



'0 r(l-a)2(l-r)2 



By a simple numerical calculation one shows that this equation does not have 
any solution for a G (0, 1). Hence, in this case there does not exist any point 
(y, a) which is an extremal of functional I[y, a]. 

b) f{y{t)) lny{t), t £ [0, 1]. Then ^ becomes 

c _ 1 

r(l-a)(l-i)" " ? 

and therefore 

y^ni^^i^tr eAci[o,i]). 

c 

In this particular case we take the set of admissible functions to be U :— {y £ 
Ui I y{l) = 0}. Using dSO]), equation ^ reads 

1 

(V'(l - a) - ln(l - t)) dr^O, 

which, after integration, yields ■(/'(a — 1) = 1. A unique solution of this equation 
in (0,1) is a = 0.604.... 

Therefore, a unique stationary point of 

I[y,a] = [ {c-oD^y + lny{t))dt 



is the point {y, a) = (no|M(i - f)"; 0, 604). 
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Remark 5.7 So far, we have considered variational problems defined via func- 
tionals of type ([T]). In fact, we have allowed fractional derivatives of functions 
to appear in Lagrangians. The natural generalization of such problems consists 
of replacing the Lebesgue integral in ([1} by the Riemann-Liouville fractional 
integral. More precisely, for /3 > set 

Ip[y.a]: = ^I^^ L{t,y{t),^Dfy,a) 



1 



(6 - v{t),oDfy, a) dt, t € (0, b). 



m Jo 

Then the fractional variational problem consists of finding extremal values of 
the functional /^[y, a]. 

In the above construction we have used the left Riemann-Liouville fractional 
integral of order /? (which, in general, differs from the order of fractional differ- 
entiation a), evaluated at t = b. The choice /3 — I turns us back to the problem 

O- 

The study of such fractional variational problems is reduced to the case 
we have already considered in the following way. It suffices to redefine the 
Lagrangian as 

L,{t, y{t),oDfy, a, (3) " tf^'W, oD^y, a). 

Then we have to consider the functional 

Ip[y,a]= [ Li{t,yit),oD^y,a,/3)dt. (32) 
Jo 

In case P > 1, Li is of the same regularity as L, so the straightforward appli- 
cation of the results derived in previous sections to the Lagrangian Li leads to 
the optimality conditions for the variational problem defined via the functional 
(|32)) . However, when < /3 < 1, continuity, as well as differentiability of Li 
with respect to t may be violated (which depends of the explicit form of L), and 
hence it may be not possible to use the theory developed so far. 

6 Conclusion 

We formulated Hamilton's principle so that the order of derivative in the La- 
grangian is also subject to variation. The stationarity conditions are derived in 
© and (ITUl) . We introduced additional assumptions which resulted in equivalent 
problems, simpler for solving. Several examples are given in order to illustrate 
the theory presented in the paper. We concluded our work with a consider- 
ation of Hamilton's principle defined in terms of Riemann-Liouville fractional 
integrals. 
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